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When  a  brittle  solid  is  loaded  to  failure,  it  docs  so  by  the  propagation  of 
cracks.  The  cracks  nucleate  and  propagate  from  inhomogeneities,  by  which  we 
mean  holes,  inclusions,  microcracks,  surface  scratches  or  other  defects.  The 
difference  between  compressive  and  tensile  fracture  is  that  in  tension  a  single 
crack  grows  unstably  (once  started,  it  accelerates  across  the  sample  to  cause 
failure)  while  in  compression  a  population  of  small  cracks  extend  stably,  each 
growing  longer  as  the  stress  is  raised,  until  they  interact  in  some  cooperative  way 
to  give  final  failure  (Figure  1).  Because  of  this,  the  strength  of  a  brittle  solid  in 
compression  is  usually  greater,  by  a  factor  of  ten  or  more,  than  that  in  tension. 

Measurements  of  the  crushing  strength  of  stone,  brick  and  of  cement  must 
have  been  of  interest  to  civil  engineers  since  pre-Roman  times.  Systematic 
measurements  of  compressive  strength  really  began  about  the  middle  of  the  last 
century  (for  its  history,  see  Jaeger  and  Cook,  1976)  but  without  much  attempt  to 
understand  what  determined  it,  or  why  brittle  materials  had  useful  strength  in 
compression  but  none  to  speak  of  in  tension.  Elucidation  of  the  mechanics  of 
brittle  tensile  fracture  has  its  roots  in  the  work  of  Griffith  (1924),  Irwin(1958) 
and  others  that  followed  (sec  Knott,  1973,  for  a  review),  which  has  led  to  the 
development  of  fracture  mechanics  as  a  branch  of  engineering  design.  The 
understanding  of  compressive  brittle  fracture  is  more  recent,  and  still 
incomplete.  A  recent  scries  of  papers  and  reviews  (Griggs  and  Handin,  1960; 
Paterson,  1978;  Hallbaucr  ct  al.,  1973;  Tapponnier  and  Brace,  1976;  Wawersik  and 
Fairhurst,  1970;  Wawersik  and  Brace,  1971;  Nemat-Nasscr  and  Horii,  1982;  Newman, 
1978;  Ashby  and  Hallam,  1986;  Sammis  and  Ashby,  1986)  have  established  that  an 
isolated  crack  in  a  large  body  grows  stably  until  its  length  becomes  comparable 
with  the  dimension  of  the  body  itself;  and  that  when  many  cracks  are  present  (as 
they  always  arc  in  natural  rocks,  in  brick,  in  concrete  and  most  ceramics)  the 
cracks  grow  stably  until  their  length  is  comparable  to  their  spacing  when  they 
interact,  an  instability  develops,  and  the  sample  fails. 
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The  problem  can  be  complicated  by  time-dependent  effects  (Anderson  and 
Grew,  1976;  Martin,  1972;  Waza  cl  al.,  1980;  Sano  ct  a!.,  1981,  Costin  and  Holcomb, 
1981;  Atkinson  and  Meredith,  1987b),  which  have  at  least  two  origins.  On  the  one 
hand  crack  growth  can  be  limited  by  a  chemical  reaction,  often  with  water.  On 
the  other,  cracking  in  compression  is  associated  with  dilation;  if  the  body  is 
saturated  with  a  fluid,  then  its  flow  into  the  dilating  region  can  introduce  a  time- 
dependent  aspect  to  fracture.  In  both  cases,  a  static  load  which  does  not 
immediately  cause  failure  may  still  do  so  if  left  in  place  for  a  sufficient  length  of 
time. 

The  understanding  of  compressive  brittle  fracture  is  still  incomplete,  but 
the  mechanisms  involved  arc  much  clearer  than  a  decade  ago.  It  seems  an 
appropriate  time  to  try  to  abstract  from  the  new  observations  and  modelling  a 
simplified  description  of  compression-cracking,  basing  it  as  far  as  possible  on  the 
physical  understanding.  The  goal  is  to  develop  a  damage  mechanics  of  brittle 
solids,  from  which  the  stress-strain  response  and  an  operational  definition  of 
failure  can  be  derived  for  a  material  with  a  given  set  of  elastic  properties  and 
given  defect  population,  under  a  given  state  of  stress.  Two  attempts  to  achieve 
this  can  be  found  in  the  open  literature;  that  of  Costin  (1983,  1985),  and  that  of 
Sammis  and  Ashby  (1986).  Central  to  the  problem  is  the  relationship  between 
stress  and  crack  extension.  Costin  (1983,  1985)  postulates  a  relationship  of 
reasonable  form,  and  develops  from  it  expressions  for  the  failure  surface  which 
(with  some  adjustable  parameters)  give  a  good  description  of  the  experimental 
data  available  at  that  time,  but  the  model  is  not  based  on  a  physical  model  for 
crack  growth.  Sammis  and  Ashby  (1986)  and  Ashby  and  Hallam  (1986)  use 
methods  of  fracture  mechanics  to  develop  a  physical  model  for  crack  extension, 
which  they  use  to  plot  stress-strain  curves  for  brittle  solids  from  which  failure 
surfaces  can  be  constructed,  but  the  complexity  of  their  model  makes  the  process 
cumbersome.  In  the  present  paper  we  attempt  to  develop  a  simpler,  model-based 
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mechanics  of  brittle  compressive  fracture,  drawing  heavily  on  the  previous 
pieces  of  work. 

2.  CRACK  INITIATION  IN  COMPRESSION 

Most  brittle  solids  contain  inhomoccncities:  small  holes  or  cracks,  particles 
which  are  poorly  bonded,  or  phases  which  have  different  moduli  or  strengths 
from  those  of  the  matrix.  Any  one  of  these  can  act  as  nuclei  for  new  cracks  when 
the  solid  is  loaded. 

The  range  of  possible  nuclei  is  wide,  but  the  spectrum  of  their 
characteristics  is  probably  bracketed  by  two  extremes:  the  spherical  hole  and  the 
sharp  inclined  crack  (Figure  2).  Both  have  been  studied  experimentally  and  both 

have  been  modeled,  the  first  by  Sammis  and  Ashby  (1986)  and  the  second  by 
Nemat-Nasscr  and  Horii  (1982)  and  Ashby  and  Hallam  (1986).  In  both  cases,  the 
criterion  for  crack  initiation,  under  axisymmetric  loading  has  the  form 

°1  =  C1  a3  "  °0 

where  ci  and  oo  arc  material  properties,  orj  is  the  axial  stress,  and  <T2=<*3  the  radial 
stress  (both  positive  when  tensile,  negative  when  compressive). 

In  the  later  development  of  this  paper  we  consider  the  growth  of  crack- 
damage  from  initial,  inclined  cracks  as  in  Figure  2a.  For  this  case  (Nemat-Nasser 
and  Horii,  1982;  Ashby  and  Hallam,  1986)  cracks  initiate  when 


where  p  is  the  coefficient  of  friction  acting  across  the  crack  faces,  Kjc  is  the 

fracture  toughness  of  the  material  through  which  the  new  crack  propagates,  and 

2a  is  the  length  of  the  original  inclined  crack.  Rocks,  typically,  show  a 
coefficient  of  friction  of  about  0.6,  in  which  case  C]  =  3.1  and  3.1  K^/ViTa  . 

Crack  initiation  from  holes  (Sammis  and  Ashby,  1986)  gives  similar  values. 
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Crack  initiation  can  be  detected  in  several  ways:  by  the  start  of  acoustic  emission, 

by  the  first  non-linearity  of  the  stress-strain  curve,  by  the  dilation  of  the  sample, 

or  by  a  sudden  increase  in  internal  friction.  None  give  very  accurate  data,  but 

they  do  allow  a  test  of  cqn.  (1).  Figure  3  shows  data  for  crack  initiation  in 

Westerly  granite  obtained  by  the  first  three  techniques  (Holcomb  and  Costin, 

1986;  Brace  et  al„  1966)  plotted  on  axes  of  oi  and  03  to  allow  comparison  with  eqn. 

(1).  The  linear  relationship  gives  a  good  description  of  the  data  with  a  slope 

between  2.7  and  3.3  (corresponding  to  p  =  0.55  to  0.64)  and  an  intercept  of  70  -  79 

1/2 

MPa  (corresponding  to  a  crack  length  2a  dose  to  1  mm  when  Kjc  =  l  MPa  m  ). 

The  theory  gives  an  adequate  description  of  the  data.  It  is  used  to  describe 
the  initiation  of  damage  in  the  diagrams  shown  later.  In  each  case,  experimental 

data  arc  fitted  to  cqn.  (1)  to  give  p  and  a  (using  published  data  for  Kic).  Results  of 


this  analysis  are  summarized  in  Table  2.  In  the  computations,  it  is  convenient  to 
normalize  the  equations  by  the  quantity  Kfc/ Vrt  a  ,  giving 


Si-c,Sj-S0 


(2) 


with 

Sj  =  fnH/Kk 
S3  =  o3  YlTa/Kic 


(3) 


Once  initiated,  the  wing  cracks  (as  we  shall  call  the  crack-like  extensions 
of  the  original  flaw)  grow  longer.  During  growth,  the  stress  intensity  Kj  at  the 
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tip  of  each  wing  crack  is  equal  to,  or  exceeds,  the  fracture  toughness  Kjc  of  the 
solid.  The  condition  for  crack  advance  is  simply 

Kl>  Kic 

The  difference  between  tension  and  compression,  as  already  mentioned,  is  that 
growth  in  compression  is  stable:  each  increment  of  crack  advance  requires  an 
increment  of  load,  at  least  until  the  cracks  start  to  interact  strongly.  We  will 
assume  that  a  steadily  increasing  load  drives  the  cracks  at  a  steady  rate,  though  in 
reality  the  inhomogcncity  of  natural  materials  may  cause  them  to  extend  in  little 
jumps.  The  problem,  then,  is  to  calculate  Kj  at  the  tip  of  the  wing  cracks. 

3.1  Crack  Growth  from  Starter  Flaws  of  a  Single  Size  :  The  2-Dimensional  Case 
Figure  4  shows  an  array  of  through-cracks,  growing  in  a  linear-elastic 
medium  under  a  triaxial  stress  field  <s\  ,  03  ,  positive  when  tensile,  negative  when 
compressive.  Consider  first  the  growth  of  a  single,  isolated  crack  from  an  initial 
inclined  flaw;  interaction  comes  later.  The  upper  inset  of  Figure  4  isolates  one 
crack:  it  is  made  up  of  an  initial  crack  of  length  2a  lying  at  an  angle  y  to  the  Xj 
direction  with  two  wings,  each  of  length  1  which  (we  will  assume)  lie  parallel  to 
X 1 .  The  stress  intensity  at  the  tips  of  the  wings  is  obtained  approximately,  but 
adequately,  in  the  following  way,  based  on  the  work  of  Nemat-Nasser  and  Horii 
(1982),  Ashby  and  Hallam  (1986),  Horii  and  Nemat-Nasser  (1985,  1986),  and 
Kemcny  and  Cook  (1986). 

The  remote  field  01,03  creates  a  shear  stress  t  and  a  normal  stress  o  across 
faces  of  the  initial  crack.  The  crack  slides  (resisted  by  the  coefficient  of  friction 
p).  wedging  open  the  mouth  of  each  wing  crack  by  5  (Figure  4).  The  wedging  can 
be  thought  of  as  caused  by  forces,  F3,  parallel  to  X3,  acting  at  the  midpoint  of  the 
crack.  The  stresses  i  and  o  arc  given  by 

ai~al  .  , 

X  = - - -  s,n  2V  (4a) 


5 


(4b) 


<T  = 


o3  +  ol 


o3-o, 


cos  2y 


2  2 

F3  is  «•'  -iply  the  component  of  the  sliding  force  acting  parallel  to  X3: 


F3  =  (t  +  )io)  2a  sinvj/ 


(5a) 


o  r 


F3  =  -(AiCV  A3O3)  a  (5b) 

where  Ai  and  A3  arc  constants  which  depend  on  \y,  to  be  determined  in  a  moment. 

The  force  F3  acting  at  the  mid  point  of  a  crack  of  length  21  creates  a  stress 
intensity  tending  to  open  the  crack  (Tada  et  al.,  1985,  page  5.1)  of 


CK0i  = 


VrTi 


This  result  gives  a  good  estimate  of  the  stress  intensity  at  the  tip  of  a  whig  craok 
when  1  is  large,  but  it  breaks  down  (becoming  infinite)  when  1  is  vanishingly 
small.  The  stress  intensity  at  the  tip  of  the  initial  inclined  crack  is  mot  in  finite, 
but  can  be  calculated  exactly  as  explained  in  the  last  section.  We  overcome  this 


problem  by  introducing  an  "effective"  crack  length  (  l  +  |ia)  giving 
(K0t=  F3 


Vtc(l  + |5a) 

We  then  choose  p  so  that  (K|)i  becomes  equal  to  that  for  the  inclined  crack  when  1 
is  zero. 

Before  doing  this,  we  note  that  there  is  another  contribution  to  K]  at  the  tip 
of  the  wing  crack.  The  remote  confining  stress  03  acts  not  just  on  the  angled 

crack  but  on  the  wing  cracks  of  length  I.  In  so  doing,  it  produces  an  additional 
contribution  to  the  stress  intensity,  tending  to  close  the  crack  when  03  is 
compressive  (Tada  et  al.,  1985,  page  5.1): 


(K[)3  =  a3Vrc7 

Summing  the  two  contributions,  with  F^  given  by  cqn  (5b),  gives: 


(7) 


vk  ( 1  +  p  a ) 


AjOj Vn  a 

Jt  Vl+  (3 


■  +  (JjVti  a 


7t  Vl+p 


h-VL 


where  L  =  1  /a.  The  cracks  extend  until  Kj  becomes  equal  to  Kjc. 

The  constants  are  found  by  ensuring  that  this  equation  reduces  to  the  exact 
result  for  crack  initiation  (L  =  0)  and  matches  the  known  results  for  very  long 
cracks  (L  »  1),  given  by  Ncmat-Nasscr  and  Horii  (1982)  and  Ashby  and  Hallam 
(1986),  eqns.  (3)  and  (6).  This  gives 


V3 


A3  -  At 


(  2\  U2 

l  1  +  M-  )  +\l 

f  x  1  /2 

(  1  +  |a2)  -\i 


p  =  0.1 

Equation  (8b)  with  the  values  of  A],  A3,  and  P  (and  with  Ki  =  K(c  )  is  plotted  in 
Figure  5.  It  shows  Ctj  Vk  a  /  K  Ic  plotted  against  L  with  the  earlier  numerical 

results  of  Ncmat-Nasscr  and  Horii  (1982).  Equation  (8)  is  obviously  a  good 
approximation  to  the  earlier  calculations. 

Now  the  interaction.  The  main  part  of  Figure  4  shows  an  array  of  Na 
cracks  per  unit  area,  all  of  which  have  extended  to  a  length  2(  1  +  a  a.)  The  center 


to-ccntcr  spacing  of  the  cracks  is 


so  that  an  uncracked  ligament  of  average  length  S  -  2  (  1  +  a  a)  remains  between 
the  cracks  in  the  Xj  direction.  (Here  a  is  simply  a  geometric  constant,  and  must  be 
distinguished  from  (5;  for  cracks  at  45°  toXi,a=l/V2).  An  opening  force  F3 

acts  at  the  midpoint  of  each  crack.  Equilibrium  requires  that  this  opening  force 
be  balanced  by  a  mean  internal  stress  03*  in  the  matrix,  as  shown  in  the  right 
hand  side  of  Figure  4.  The  average  internal  stress  is  given  by 


a3  = 


This  acts  on  the  wing  cracks 


S  -2(l  +  aa) 


(11) 


(K03  =  (o3+a3)vTl  (12) 

We  now  define  the  initial  damage  Dq  and  the  current  damage  D  by: 


D0  =  Jt  (  a  a  )  Na 


(13a) 


D  =  n  ( 1  +  a  a)  NA 


giving 


i  -(AjtJj-  A3a3)(D0/rc) 


1/2 


<*3  = 


a  (l  -2(D/K)m) 

Equations  (8a)  and  (8b)  now  become 

Ki  = - — - +  (a,  +  a,  )VtTi 

Oca+ea))1'2 


(13b) 


(14) 


(15a) 
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Here  the  first  term  in  the  curly  brackets  describes  the  wedging  plus  the  crack- 
crack  interaction;  and  the  second  term  describes  the  closing  effect  of  the  lateral 
confining  stress.  The  cracks  propagate  until  Ki  falls  to  Kic-  Using  this, 

rearranging  and  aggregating  the  constants  (with  y  2/ K  ~  1  )  gives  for 

proportional  loading  (with  X  =  03/01  held  constant)  : 


and  for  loading  at  constant  03 


(17) 

where  Sj  and  S3  arc  defined  by  equation  (3).  The  values  of  the  constants  are 
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(18) 


Figure  6  shows  how  the  axial  stress  <T]  varies  with  damage  D  for  various  confining 
pressures.  The  left  hand  figure  shows  proportional  loading;  the  right  hand 
figure  ,  loading  at  constant  03.  The  peak  stress,  (oi)max  .  rises  and  moves  to  the 
right  as  X  or  03  is  increased.  The  shapes  of  the  curves  at  constant  X  differ  from 
those  at  constant  03,  as  expected,  but  the  peaks  are  at  the  same  stress.  Figure  7 
shows  (<ri)max  plotted  against  03  for  both  conditions;  the  points  lie  on  the  same 
line. 

3.2  Crack  Growth  from  Starter  Flaws  of  a  Single  Size.  The  3-dimcnsional  Case 
It  is  usually  the  case  that  flaws  are  completely  contained  within  the 
material.  The  merit  of  the  2-dimcnsional  calculations  developed  in  section  3.2  is 
that  it  points  to  a  way  of  tackling  this  more  difficult  3-dimensional  problem.  We 
require  the  stress  intensity  at  the  periphery  of  a  contained  crack  emanating  from 
a  starter  flaw  (which  we  take  to  be  an  inclined,  penny-shaped  crack)  which  we 
will  equate,  as  before,  to  Kjc  .  This  we  do  by  calculating  the  wedging  force  F3  as 
in  the  2-dimensional  case.  The  wedging  force  creates  an  average  internal  stress 
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03 1  •  The  stress  intensity  at  the  tip  of  a  given  wing  crack  is  calculated  from  the 
wedging  force  and  the  total  lateral  stress  (03  +03'),  as  shown  in  Figure  8.  The 
significant  difference  in  the  results  of  this  3-D  and  the  earlier  2-D  calculations  is 
that  the  dependence  of  Si  and  S3  on  damage  D  involves  different  powers. 

The  wedging  force  F3,  as  before,  is  calculated  from  the  shear  and  normal 
stresses  (cqn  4)  acting  on  the  initial  crack  plane,  times  the  crack  area,  resolved 
into  the  X3  direction: 

2 

F3  =  (t+4  0)7ta  siny 

=  -  (A,  o,  -  A3 ct3  )  a2  (19) 

Unlike  the  2-D  case,  there  arc  no  exact  analytical  solutions  for  limiting  cases 
which  allow  Ai  and  A3  to  be  determined,  so  we  make  the  assumption  that  they 
have  the  same  values  as  before  (equation  9)  but  make  provision  to  adjust  them 
later  to  match  experimental  data. 

The  stress  intensity  Ki  has  the  same  three  contributions  as  before.  The 
wedging  force  F3  induces  a  stress  intensity  (Tada  ct  al.,  1985,  page  24.2) 

f3 

,  x3/2  (20) 

[n  (  1  +  p  a  )J 

where  P  is  introduced  for  the  same  reason  as  before:  to  give  a  limiting  value  of 
( K 1)  1  when  I  =  0.  The  contribution  due  to  a 3  and  03'  arc  (Tada  et  al.,  1985,  page 
24.2). 


The  internal  stress  for  the  3-D  case  is 


°3  = 


A  - n(l  +  a  a) 


1  1 


(21) 


(22) 


where  JC  (l  +  a  a)  is  the  total  crack  area  projected  normal  to  X3  and  A  is  the 
per  crack 


area 


A  =  71 


"3/  3 


2/3 


4  N 1 


(23) 


where  Ny  is  the  number  of  cracks  per  unit  volume. 

Damage  is  defined  in  a  way  which  parallels  that  in  2-dimcnsions: 


4  3 

D0  =  y  71  ( a  a )  N, 


(24a) 


4  3 

D  =  y7C(l  +  aa)  N, 


giving 


i  -  (  A  j  <Tj  -  A3  Oj  )  Dq 


2/3 


o3  = 


7t  a  (l  -  D2/3) 

The  stress  intensity  at  the  tip  of  the  wing  crack  is 


K,  = 


r  7i 

(tc  ( I  +  P  a ) ) 


-  (a3  +  a3)  VtTi 


(24b) 


(25) 


(26) 


-  A,  a,  Vn~a  ff  A,  \ 

n2a3'2((D»0)''3.l+p/«)5'2^  A'  ] 


2/  „  W 


1  + 


2((EVD0),/3-l)  -5® 


1  -D 


2/3, 


•77a2*2((D/D«>)''5-1)2} 


As  before  the  cracks  propagate  until  K{  fails  to  Kjc.  Rearranging  and 
aggregating  the  constants  the  gives,  for  proportional  loading: 


The  equations  and  constants  have  a  form  very  like  those  of  the  2- 
dimensional  model.  Two  significant  differences  should  be  noted.  First,  the  extra 
dimension  causes  the  powers  of  D  which  appear  in  the  equation  to  differ  (not 

surprisingly)  from  those  of  the  2-D  model.  Second,  the  constants  Cj  to  C4  are  not 

known  with  the  same  precision  as  (hose  of  the  2-D  model  because  accurately- 
known  limiting  cases  arc  not  available  to  calibrate  them.  We  shall  assume 
(reasonably)  that  the  dependence  on  the  coefficient,  p,  is  properly  included,  but 

that  the  constant  p  may  require  further  adjustment  to  give  a  good  match  with 
experiment. 

Figures  9  and  10  illustrate  some  features  of  the  results.  The  axial  stress  at 

first  rises  as  damage  grows  (Figure  9),  passing  through  a  peak  which  shifts  to 

higher  values  of  damage  as  the  confining  pressure  increases.  The  damage 

surface,  shown  in  Figure  10,  is  almost  a  cone  meaning  that,  to  a  first 

approximation,  the  failure  envelope  is  described  by 

o,  =  Co3  -  ac 

where  C  is  a  constant  and  oc  is  the  unconfincd  compressive  strength.  The  value  of 
the  model  is  that  it  gives  a  physical  interpretation  to  C  and  ac  ,  and  relates  them  to 
the  initial  damage,  the  coefficient  of  friction,  and  the  crack  size. 

4.  ANALYSIS  OF  DATA 

The  strength  of  many  different  rocks  have  been  determined  under  Biaxial 
loading  conditions  (  la ]  1  >  Ia2l  =  l<T3 1,  where  all  three  principal  stresses  am 
compressive)  .  For  a  few,  the  initiation  of  microcracking  has  also  been 
determined.  We  now  apply  the  damage  mechanics  model  developed  above  to  nine 
different  rock  types  for  which  the  most  complete  data  sets  exist  :  granite,  tplite, 
dunitc,  eclogitc,  gabbro,  sandstone,  limestone,  marble  and  rock  salt.  These  rocks 

represent  a  wide  range  of  composition  and  initial  damage.  Granite,  aplifev  dunitc 
and  cclogite  arc  low  porosity,  crystalline,  igneous  rocks  in  which  the  iaitiaf 
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damage  is  mostly  in  the  form  of  low  aspect  cracks.  Limestone  and  sandstone  arc 
porous  sedimentary  rocks  in  which  the  initial  damage  is  mostly  in  the  form  of 
high  aspect  pores.  Marble  is  a  metamorphic  rocks  with  initial  damage  of  a  form 
intermediate  between  the  previous  two  extremes.  These  rocks  also  span  a  wide 
range  of  yield  strength.  The  igneous  rocks  have  yield  strengths  in  excess  of  2  GPa 
while  the  calcareous  rocks  (limestone  and  marble)  yield  at  stresses  below  1  GPa. 
Rock  salt  is  at  the  low  extreme  with  a  yield  strength  below  100  MPa. 

The  damage  mechanics  model  formulated  above  has  three  constants:  p  a, 
and  p.  They  arc  not  strictly  adjustable  because  the  theory  makes  predictions  for 
their  values.  However  because  of  approximations  in  the  derivation,  and  the 
uncertainties  in  the  aspect  ratio  of  the  starter  flaws,  we  have  treated  p  as 

adjustable,  choosing  the  value  0.45  to  give  the  best  fit  to  the  data. 

For  those  materials  where  crack  initiation  data  are  available,  the  crack 
length  2a  and  the  coefficient  of  friction  p  are  determined  from  the  initiation 
surface  (see  Table  2) 

<j\  =  Ci  03  +  o0 

where  Cl  and  o0  are  given  by  equation  (1).  Kic  is  also  required.  Although  Kjc 

may  be  estimated  for  most  rocks  (Atkinson  and  Meredith,  1987a),  the  starter  flaw 
size  is  not  known  in  most  cases  and  must  be  treated  as  an  adjustable  parameter. 

The  derivation  of  the  fundamental  equations  (27)  and  (28)  from  equation  (26) 
gives  C3  =  2,  and  this  gives  a  good  description  of  the  materials  we  have  examined. 

Each  material  will  now  be  discussed  in  tum.  The  triaxial  data  for  damage 
initiation  and  failure  arc  presented  on  plots  of  01  vs  03  .  The  theoretical  fracture 
initiation  surface  (cqn.  1),  surfaces  of  constant  damage,  and  the  failure  surface 
(calculated  from  the  maximum  of  cqn.  (28)  for  each  value  of  <73  )  are  plotted  on 
these  graphs  for  comparison  with  the  data. 

As  the  confining  pressure  is  increased,  brittle  fracture  is  made 
increasingly  difficult.  A  critical  pressure  may  be  reached  at  which  true  plasticity 


replaces  crack  extension.  This  transition  can  be  illustrated  by  plotting  a  yield  (or 
creep)  surface,  defined  by: 


2  i  I"  2  2  2*1 

ay  =  j  |(<Vct2>  +(<Va3)  +  (<T3  -  <t,  )  J  (30) 

The  yield  surface  is  plotted  as  a  pair  of  heavy  broken  lines  on  each  figure  The 
yield  strength,  Oy,  can  be  derived  from  hardness,  H,  data  since  <Ty  =  H/3.  The 
material  properties  and  constants  used  to  generate  the  theoretical  initiation  and 
failure  curves  arc  tabulated  for  each  solid. 

4.1  Qianiic 

Westerly  granite  is  a  fine-grained  (0.75  mm.),  low  porosity  (0.9%), 
isotropic,  two-mica  calc-alkaline  granite  which  has  become  a  standard  material  in 
rock  mechanics  testing  (sec  Scholz,  1986.  for  a  brief  history).  Mineralogical 
modal  analyses  arc  given  by  Birch  (1960)  and  Wawersik  and  Brace  (1971). 

Figure  11  a  shows  theoretical  surfaces  for  initiation,  constant  damage,  and 
failure  at  low  values  of  the  confining  stress.  The  fracture  initiation  data  from 
Brace  et  al.  (1966)  were  determined  from  the  onset  of  non-linear  behavior  of  the 
volume  strain.  Only  data  taken  at  the  highest  loading  rate  are  plotted  here  in 
order  to  minimize  effects  of  subcritical  crack  growth  (which  we  do  not  model). 

The  three  initiation  points  from  Holcomb  and  Costin  (1986)  were  determined  from 
the  onset  of  acoustic  emission  (AE)  in  a  previously  unstressed  sample.  Also  shown 
is  a  surface  of  constant  damage  mapped  by  Holcomb  and  Costin  (1986)  using  the  AE 
Kaiser  effect  as  a  probe.  The  triaxial  failure  data  are  from  Brace  et  al.  (1966)  and 
Mogi  (1966). 

Figure  lib  shows  the  failure  surface  extended  out  to  large  values  of  the 
confining  stress  where  it  intersects  the  yield  surface.  Data  at  low  and 
intermediate  confining  pressures  are  from  the  same  sources  as  in  Figure  11a. 
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Those  at  high  confining  pressures  are  from  Schock  and  Heard  (1974)  for  Westerly 
granite  and  from  Shimada  (1981)  for  Man*nari  granite  (grain  size  1-3  mm. 
apparent  porosity:  0.7%).  It  is  evident  that  the  failure  surface  has  considerable 
curvature  and  deviates  from  our  theoretical  model  at  high  confining  pressures. 
Although  Janach  and  Guex  (1980)  have  modeled  this  curvature  in  terms  of  the 
formation  of  shear  bubbles  at  the  grain  boundaries.  Figure  lib  supports  the 
possibility  that  the  curvature  is  due  to  a  gradual  transition  to  ductile  behavior. 
Analogous  curvature  is  evident  in  subsequent  figures  for  limestone,  marble,  and 
NaCl  which  arc  known  to  exhibit  ductile  behavior  at  moderate  confining 
pressure,  although  the  curvature  in  these  rocks  occurs  over  a  more  limited 
pressure  range.  The  broader  transition  in  granite  may  reflect  its  multimineralic 
composition  for  which  the  individual  minerals  have  different  brittle  -  ductile 
transition  pressures.  Note  that  the  two  granites  show  different  transitional 
behavior  at  intermediate  pressures,  but  that  both  approach  the  same  ductile  limit. 
Schock  and  Heard’s  observation  that  the  stress-  strain  curve  is  linear  to  failure 
may  be  due  to  the  convergence  of  the  initiation  and  failure  surfaces  at  very  high 
confining  pressures  in  Figure  lib.  The  sudden  release  of  energy  and  shear 
localization  at  failure  do  not  preclude  stress  concentration  by  ductile  processes  in 
the  weaker  minerals. 

4.2  Aplite 

Brace  ct  al.  (1966)  studied  a  quartz-oligoclasc  aplite  (63%  oligoclase,  27% 
quartz,  and  10%  biotite)  which  they  described  as  fine-grained  and  flinty, 
apparently  isotropic,  and  of  high  strength.  The  feldspar  is  highly  altered.  The 
grain  size  of  the  ground  mass  is  about  40  pm  and  of  the  phcnocrysts,  about  100  pm. 
The  small  grain  size  and  high  strength  arc  consistent  with  the  small  flaw  size 
required  by  the  model  (sec  Figure  12).  The  relatively  high  initial  damage  is 
consistent  with  the  flinty  texture.  A  high  density  of  small  flaws  may  explain  why 


flinty  materials  can  be  reliably  fashioned  into  tools  by  flaking  off  small  bits  in  a 
controllable  manner. 

4.3  Dunite 

Dunite  is  an  almost  pure  olivine  rock.  Shimada  et  al.  (1983)  measured  the 
compressive  strength  of  Horoman  dunite  (grain  size  0.1 -0.9  mm)  at  confining 
pressures  up  to  450  MPa  using  a  conventional  triaxial  testing  apparatus  and  to  3 
GPa  using  a  cubic  press.  Acoustic  emissions  showed  a  change  in  failure  mode  at 
confining  pressures  between  0.44  and  0.77  GPa.  Below  these  pressures,  AE  activity 
increased  rapidly  between  the  onset  of  dilatancy  and  failure  (the  typical  pattern 
for  brittle  failure).  Above  the  transition  pressure,  and  increase  in  AE  activity  was 
not  observed  to  precede  failure;  rather,  the  level  of  AE  remained  nearly  constant 
up  to  failure.  Shimada  cl  al.  (1983)  correlate  this  change  in  AE  behavior  with  the 
extreme  curvature  in  the  failure  envelope.  As  is  evident  in  Figure  13,  our  model 
suggest  that  this  change  in  behavior  is  associated  with  the  transition  to  plastic 
deformation. 

4.4  Eclogitc 

Eclogite  is  an  ultramafic  pyroxenr-gamet  rock.  The  Akaishi  eclogite 
measured  by  Shimada  et  al.  (1983)  was  composed  of  0.1 -0.3  mm  pyroxene  grains 
and  0.8-2. 3  mm  garnet  grains.  It  had  a  density  of  3.642  gm/cm^  and  a  porosity  of 
0.4%.  Conventional  triaxial  tests  covered  a  range  of  confining  pressures  from  0  to 
450  MPa  while  tests  in  an  opposed  anvil  cubic  press  extended  the  confining 
pressure  to  3  GPa.  As  discussed  above  for  dunite,  the  AE  patterns  indicate  a 
change  in  failure  mechanism  at  confining  pressures  between  1.02  and  1.99  GPa. 
The  data  and  theoretical  surfaces  arc  given  in  Figure  14. 

4.5  Gabbro 

In  addition  to  dunite  and  eclogite  discussed  above,  Shimada  et  al.  (1983)  also 
studied  Murotomisaki  gabbro,  a  hypersthcnc-bcaring-olivinc-augite  gabbro.  The 
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grain  size  of  the  olivine  component  is  1-2  mm,  pyroxene  is  about  0.7  mm,  and  the 
plagioclasc  is  about  0.7-3  mm.  The  bulk  density  is  2.985  gm/cm^  and  the  reported 
porosity  is  0.4%.  The  pattern  of  AE  indicated  a  change  in  failure  mechanism  at 
confining  pressures  between  0.51  and  0.76  GPa.  The  data  and  theoretical  surfaces 
are  given  in  Figure  15. 

4.6  Sandstone 

The  sandstone  data  in  Figure  16  were  obtained  in  triaxial  compression  by 
Gowd  and  Rummcl  (1980).  The  rock  is  described  as  a  medium  grain-sized 
Buntsandstonc  from  SW-Germany  with  subangular  to  round  quartz  grains  bedded 
within  a  clayey  matrix..  Us  initial  porosity  was  15%  with  an  initial  permeability 

of  50  microdarcy.  The  damage  initiation  data  were  defined  by  the  onset  of 
dilatancy.  At  confining  pressures  above  about  30  MPa,  the  stress-strain  curves 
are  nonlinear  at  lower  values  of  the  axial  stress  than  the  observed  onset  of 
dilatancy.  This  probably  reflects  the  suppression  of  dilatancy  by  pore  collapse,  a 

phenomenon  which  the  authors  propose  to  explain  the  total  lack  of  observed 
dilatancy  at  the  highest  confining  pressures.  Such  effects  are  beyond  the  scope  of 
our  model. 

A  transition  from  brittle  failure  to  apparent  ductile  shear  deformation 
takes  place  at  a  pressure  of  about  100  MPa.  However,  the  observed  pressure 
dependence  of  the  flow  stress  for  confining  pressures  in  excess  of  100  MPa  argues 
against  true  ductile  flow  and  for  a  cataclastic  mode  of  deformation  probably 
involving  pore  collapse.  Dilatancy  at  failure  is  a  constant  for  confining 

pressures  between  0  and  40  MPa.  From  40  to  100  MPa,  dilatancy  at  failure  decreases 
to  zero.  Above  100  MPa,  brittle  failure  docs  not  occur. 

4.7  Limestone 

Solcnhofen  limestone  is  a  fine  grained  (0.01  mm)  mechanically  isotropic 
limestone  from  Bavaria.  It  has  a  connected  porosity  of  5.3%  and  a  total  porosity  in 
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the  range  6  -  9%  (Rutter,  1972).  The  strength  data  in  Figure  17  are  from  Heard 
(1960)  and  include  both  triaxial  compression  and  tension.  The  fracture  initiation 
points  were  picked  as  the  onset  of  nonlinearity  in  his  published  stress-strain 
curves  and  are  only  approximate. 

4.8  Marble 

The  only  marble  for  which  fracture  initiation  data  is  available  is  described 
by  Brace  et  al.  (1966)  as  a  medium  grained  almost  pure  calcitc  marble  of  unknown 
origin.  They  report  it  to  be  apparently  isotropic,  very  ductile  even  at  low 
confining  pressures,  and  having  a  grain  si/c  of  about  0.2  mm. 

The  unusually  low  fracture  initiation  stress  (Figure  18)  requires  either 
large  starter  flaws  (~  6  mm)  or  a  low  fracture  toughness.  Since  Atkinson  and 
Meredith  (1987a)  report  Kic  as  low  as  0.19  MPa  m^  for  calcite,  we  have  fit  the 
initiation  data  using  this  value  which  then  implies  a  starter  flaw  size  of  0.4  mm, 
which  is  comparable  to  the  grain  size. 

In  a  microscopic  study  of  nucleation  in  marble,  Olsson  and  Peng  (1976) 
found  that  microcracks  often  nucleate  where  slip  bands  intersect  grain 
boundaries.  Although  such  slip-bands  arc  physically  analogous  to  angle  cracks, 
there  may  be  a  significantly  larger  number  of  such  nuclei  since  every  favorably 
aligned  grain  is  a  potential  source  of  nuclei.  This  may  explain  the  large  values  of 
initial  damage  Do  required  to  fit  the  marble  data. 

The  data  set  which  we  fit  is  for  Carrara  marble  which  is  the  fine  grained 
(about  0.1  mm)  isotropic  marble  used  by  Michaclangclo  for  the  Picta  and  other 
well-known  works.  Its  total  porosity  is  about  1.1%  (Edmond  and  Paterson,  1972). 
The  triaxial  data  in  Figure  18  arc  from  Von  Karman  (1911)  and  Edmond  and 
Paterson  (1972).  Brittle  v&  ductile  behavior  was  deduced  from  the  shape  of  the 
stress-strain  curve  and  from  the  volume  changes  associated  with  the 
deformation. 
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Rock  salt  exhibits  a  room-temperature  brittle-to-ductilc  transition  at  the 
lowest  confining  pressure  of  any  rock  in  this  study.  Hunsche  (1981)  tested  three 
types  of  natural  salt  at  low  confining  pressure  under  both  the  common  triaxial 
loading  and  the  less  common  true  multiaxial  loading  at  strain  rates  of  about  10*6  s- 
Handin  (1953)  collected  conventional  triaxial  data  to  higher  confining 
pressures  at  a  strain  rate  of  about  10*^  s'*.  These  data  are  shown  in  Figure  19. 

5.  CONCLUSIONS 

1.  An  approximate  physical  model  for  damage  evolution  in  brittle  solids 
under  compressive  stress  states  has  been  developed.  The  model  is  based  on  the 
growth  of  wing  cracks  from  a  population  of  small,  inclined,  starter  cracks;  and 
the  interaction  between  them.  The  important  variables  of  the  problem  are:  the 
size,  2a,  of  the  initial  inclined  cracks,  and  the  initial  damage 

4  3 

D  =  jTta  Nv 

The  state  of  the  material  is  measured  by  the  current  value  of  the  damage 

D  =  ^7i(l  +  aa)  Nv 

where  1  is  the  length  of  the  wing  cracks. 

A  number  of  simplifying  assumptions  arc  made  in  the  current  model.  In 
particular  it  is  assumed  that  the  population  of  initial  cracks  all  have  the  same  size. 
One  consequence  of  this  is  that  the  initiation  surface  (that  is  ,  the  combination  of 
stresses  corresponding  to  the  first  extension  of  the  cracks,  and  thus  the  first 
increment  of  new  damage)  is  linear,  described  by: 

o':  =C]  03  -o0 


Similarly,  damage  itself  is  a  state  variable  so  that  surfaces  of  constant  damage  in 
stress  space  arc  also  linear.  The  surfaces  corresponding  to  final  macroscopic 
fracture  is  not  one  of  constant  damage  (as  often  assumed)  .  The.  terminal  damage 


itself  depends  on  the  stress  slate,  but  it,  too,  is  well  approximated  by  a  linear 
relationship: 

°1  =  C2  <*3  •  °c 

where  oc  is  the  simple  compressive  strength. 

The  model  has  been  fitted  to  data  for  a  number  of  rocks.  The  process  gives 
physical  insight  into  the  damage  accumulation  and  failure  of  these  materials  in 
compression.  In  particular,  the  fitting  process  leads  to  a  value  for  the  coefficient 
of  friction  across  the  crack  faces;  the  size  of  the  initial  flaws;  and  the  initial 
damage  Dq  .  The  failure  process  depends  principally  on  these  variables. 

Curvature  of  the  failure  surface  is  shown  to  depend,  at  least  partly,  on  an 
interaction  between  the  brittle  failure  mechanism,  and  plastic  flow.  Rocks  which 
show  clearly  established  plasticity  at  high  pressures  (marbles,  and  rock  salt,  for 
example)  show  a  brittle  regime  at  low  pressures,  a  transitional  regime  at 
intermediate  pressures  (both  depending  strongly  on  pressure),  and  a  regime  of 
plasticity  at  high  pressures  which  is  independent  of  pressure  itself.  It  is 

noteworthy  that  silicate  rocks  such  as  granite,  gabbro,  dunitc  and  aplite  show  a 
similar  behavior,  with  the  transition  to  plasticity  dominating  failure  at  confining 

pressures  of  general  order  E/30.  This  transition,  at  first  sight  a  surprising  one,  is 
nonetheless  to  be  expected  at  such  stress  levels  which  are  roughly  the  theoretical 

shear  strength  of  the  minerals  within  the  rock. 

Several  other  noteworthy  conclusions  emerge.  One  is  that,  in  rocks  which 
are  almost  fully  dense,  the  initial  flaw  size  is  roughly  equal  to  the  grain  size  of  the 
rock  itself.  But  the  initial  damage  level,  Dq  ,  varies  widely.  In  low  porosity 
crystalline  silicates  such  as  granite,  this  level  is  low  (typically  3%);  but  in 
intrinsically-plastic  materials  like  calcitc  and  rock  salt,  the  initial  damage  level  is 
high  (of  order  15%)  perhaps  because  the  flaws  from  which  wing  cracks  grow  arc 
slip  bands  within  suitably  oriented  grains  rather  than  cracks. 

3.  Data  for  rocks  which  are  almost  fully  dense  arc  well  fitted  by  the  model. 

We  find,  too,  that  porous  rocks  (limestone  and  sandstone,  both  with  roughly  15% 
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porosity)  are  also  well  fitted.  This  suggests  that  an  analogous  theoretical 
development  may  be  possible  for  porosity  induced  cracking  too. 
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Figure  1.  Compressive  failure  of  a  brittle  solid  containing  a  distribution  of 

flaws. 

Figure  2.  Cracks  can  initiate  at  inclined  flaws  and  at  holes.  In  both  cases  there 
are  two  contributions  to  K[  ,  the  opening  stress  intensity  at  the  tip  of 
the  growing  wing  cracks.  One  is  caused  by  the  stress  concentrations 
at  the  flaw;  the  other  is  due  to  the  closing  effect  of  03. 

Figure  3.  Data  for  crack  initiation  in  Westerly  granite.  Crack  initiation  data 

for  several  other  rocks  are  analyzed  in  section  4.  In  all  cases  the 
data  are  well  fitted  by  equation  (1)  with  a  coefficient  of  friction 
between  0.55  and  0.65. 

Figure  4.  A  population  of  growing  cracks.  We  first  analyze  the  growth  of  an 

isolated  crack  (shown  above)  and  then  include  the  crack-crack 
interaction  (illustrated  on  the  right). 

Figure  5.  A  comparison  of  the  approximate  equation  (8b)  with  the  numerical 

calculations  of  Ncmat-Nasser  and  Horii  (1982).  The  approximation  is 
adequate  for  the  present  purposes. 

Figure  6.  The  dependence  of  axial  stress  on  damage  as  predicted  by  equations 

(16)  and  (17).  The  peak  stress  is  marked.  We  take  this  as  the  failure 
stress. 

Figure  7.  The  peak  value  of  «i  plotted  against  03  to  show  the  failure  surface. 

In  the  present  formulation,  damage  is  a  state  variable. 

Figure  8.  Wing  cracks  growing  from  an  initial,  constrained,  penny-shaped 

flaw.  The  geometry  is  more  complicated  than  in  the  2-dimensional 
case  but  the  same  method  can  be  used  to  give  an  approximate 
solution  for  Ki  at  the  tips  of  the  wing  cracks. 

Figure  9.  The  dependence  of  axial  stress  on  damage  as  predicted  by  equations 

(27)  and  (28).  The  peak  stress  is  marked.  We  take  this  as  the  failure 
stress. 

Figure  10.  The  peak  value  of  ai  plotted  against  03  to  show  the  failure  surface. 

In  this  formulation,  damage  is  a  state  variable. 

Figure  11.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  granite.  Data  and  theory  for  microfracture  initiation  and 
surfaces  of  constant  damage  arc  also  compared.  The  yield  surface 
(cqn.  30)  is  also  plotted  as  the  heavy  broken  lines.  Part  (a)  is  limited 
to  low  and  intermediate  confining  stress;  part  (b)  extends  to  the 
largest  measured  confining  stress. 

Figure  12.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  aplitc.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 


constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 
30).  Part  (a)  is  limited  to  low  and  intermediate  confining  stress;  part 
(b)  extends  to  the  largest  measured  confining  stress. 

Figure  13.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  dunite.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30). 

Figure  14.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  eclogitc.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30). 

Figure  15.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  gabbro.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30). 

Figure  16.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  sandstone.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 
30). 

Figure  17.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  limestone.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  fight  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30). 

Figure  18.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  marble.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30).  Parts  (a)  and  (b)  are  limited  to  low  and  intermediate  confining 
stress;  part  (c)  extends  to  the  largest  measured  confining  stress. 

Figure  19.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  rock  salt.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  arc  surfaces  of 

constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30).  Part  (a)  is  limited  to  low  and  intermediate  confining  stress;  part 
(b)  extends  to  the  largest  measured  confining  stress. 
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TAB!E  1.  MATERIAL  PROPERTIES  OF  THE  ROCKS  AND  FITTED  MODEL  PARAMETERS  (Bracketted  data  are  estimates). 
MATERIALS)  K  (GPa)  £  (GPa)  Sy  (MPa)  Kic  (MPam  )  2a  (mm)  Mu  Do  C3 
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Compressive  failure  of  a  brittle  solid  containing  a  distribution  of  flaws. 


t  t  ",  t  t  t  u  t  t 


Figure  2.  Cracks  can  initiate  at  inclined  flaws  and  at  holes.  In  both  cases  there  are 
two  contributions  to  K| ,  the  opening  stress  intensity  at  the  tip  of  the 
growing  wing  cracks.  One  is  caused  by  the  stress  concentrations  at  the 
flaw;  the  other  is  due  to  the  closing  effect  of  03. 
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RADIAL  STRESS  cr3  (GPa) 


Data  for  crack  initiation  in  Westerly  granite.  Crack  initiation  data  for 
several  other  rocks  are  analyzed  in  section  4.  In  all  cases  the  data  are  well 
fitted  by  equation  (1)  with  a  coefficient  of  friction  between  0.55  and  0.65. 


NORMALISED  STRESS  -  cr1-/7ta/Klc  NORMALISED  STRESS 


LOG  (D) 


Figure  6. 


The  dependence  of  axial  stress  on  damage  as  predicted  by  equations  (16) 
and  (17).  The  peak  stress  is  marked.  We  take  this  as  the  failure  stress. 


Figure  7. 


I 


I 
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The  peak  value  of  <s\  plotted  against  03  to  show  the  failure  surface.  In  the 
present  formulation,  damage  is  a  state  variable. 
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Wing  creeks  growing  from  an  initial,  constrained,  penny-shaped  flaw.  The 
geometry  is  more  complicated  than  in  the  2-dimensional  case  but  the  same 
method  can  be  used  to  give  an  approximate  solution  for  K|  at  the  tips  of  the 
wing  cracks. 


Figure  8. 
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Figure  11.  Comparison  beiween  experimental  and  theoretical  failure  surfaces 
for  granite  at  low  and  intermediate  confining  stress.  Data  and 
theory  for  microfracture  initiation  and  surfaces  of  constant  damage 
are  also  compared.  The  yield  surface  (eqn.  30)  is  also  plotted  as  the 
heavy  broken  lines. 


RADIAL 


NORMALIZED  RADIAL  STRESS  cr3/Kx10‘ 


AXIAL  STRESS  <r,  (MPa) 


Figure  12.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  granite  extended  to  the  largest  measured  confining  stress.  Data 
and  theory  for  microfracture  initiation  and  surfaces  of  constant 
damage  are  also  compared.  The  yield  surface  (eqn.  30)  is  also  plotted 
as  the  heavy  broken  lines. 
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RADIAL  STRESS  cr3  (MPa) 


NORMALIZED  RA 


AXIAL  STRESS,  a,  (MPa) 


NORMALIZED  AXIAL  STRESS,  a;  /K  xIO'3 


Figure  13.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  aplite  at  low  and  intermediate  confining  stress.  The  heavy  solid 
line  is  the  theoretical  failure  surface.  The  light  solid  line  is  the 
surface  for  the  initiation  of  microfracturing  while  the  light  broken 
lines  are  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 
yield  surface  (eqn.  30). 
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NORMALIZED  RADIAL  STRESS,  cr, /Kx  10 
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AXIAL  STRESS,  O',  (MPa) 
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Figure  14.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  aplite  extended  to  the  largest  measured  confining  stress.  The 
heavy  solid  line  is  the  theoretical  failure  surface.  The  light  solid 
line  is  the  surface  for  the  initiation  of  microfracturing  while  the 
light  broken  lines  are  surfaces  of  constant  damage.  The  heavy 
broken  line  is  the  yield  surface  (eqn.  30). 
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RADIAL  STRESS,  a3  (MPa) 


NORMALIZED  RADIAL  STRESS.  ct3  /KxIO 


AXIAL  STRESS,  0}  (MPa) 


Figure  IS.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  dunite.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eon. 
30). 
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RADIAL  STRESS,  c r  (MPa) 


NORMALIZED  RADIAL  STRESS,  cr3/Kxl6 
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Figure  16.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  eclogite.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 
30). 
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Figure  17.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  gabbro.  The  heavy  solid  line  is  the  theoretical  failure  surface. 
The  light  solid  line  is  the  surface  for  the  initiation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 
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Figure  18 


AXIAL  STRESS  a,  (MPa) 
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The  light  solid  line  is  the  surface  for  the  initiation  of 
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RADIAL  STRESS  oc  (MPa) 
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Fieure  19  Comparison  between  experimental  and  theoretical  failure  surfaces 
F  gure  r=ne.  The  heavy  solid  line  is  the  theoretical  failure  surface. 

The  light  solid  line  is  the  surface  for  the  mutation  of 
microfracturing  while  the  light  broken  lines  are  surfaces  of 
constant  damage.  The  heavy  broken  line  is  the  yield  surface  (eqn. 

30). 


RADIAL  STRESS,  cr2  (MPa) 


NORMALIZED  RADIAL  STRESS,  cr3/Kx10 


( 


AXIAL  STRESS,  a\  (MPa) 


NORMALIZED  AXIAL  STRESS,  C^/KxlO'3 


«- »■  s"  foirX^,  y 

fines'* are  surfaces  of  cons.an.  damage.  The  heavy  broken  toe  »  .he 
yield  surface  (eqn.  30). 
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RADIAL  STRESS,  cr3  (M Pa) 
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Figure  21.  Comparison  between  experimental  and  theoretical  failure  surfaces 

for  marble  at  low  and  intermediate  confining  stress.  The  heavy  solid 
line  is  the  theoretical  failure  surface.  The  light  solid  line  is  the 
surface  for  the  initiation  of  microfracturing  while  the  light  broken 
lines  are  surfaces  of  constant  damage.  The  heavy  broken  line  is  the 
yield  surface  (eqn.  30). 
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Figure  22.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  marble  extended  to  the  highest  measured  confining  stress.  The 
heavy  solid  line  is  the  theoretical  failure  surface.  The  light  solid 
line  is  the  surface  for  the  initiation  of  microfracturing  while  the 
light  broken  lines  are  surfaces  of  constant  damage.  The  heavy 
broken  line  is  the  yield  surface  (eqn.  30). 
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RADIAL  STRESS,  or,  (MPa) 


NORMALIZED  RADIAL  STRESS.  cr3/ Kx10~ 


Figure  23.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  rock  salt  at  low  and  intermediate  confining  stress.  The  heavy 
solid  line  is  the  theoretical  failure  surface.  The  light  solid  line  is 
the  surface  for  the  initiation  of  microfracturing  while  the  light 
broken  lines  are  surfaces  of  constant  damage.  The  heavy  broken 
line  is  the  yield  surface  (eqn.  30). 
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RADIAL  STRESS.  cr3  (MPa) 
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Figure  24.  Comparison  between  experimental  and  theoretical  failure  surfaces 
for  rock  salt  extended  to  the  largest  confining  stress.  The  heavy 
solid  line  is  the  theoretical  failure  surface.  The  light  solid  line  is 
the  surface  for  the  initiation  of  microfracturing  while  the  light 
broken  lines  are  surfaces  of  constant  damage.  The  heavy  broken 
line  is  the  yield  surface  (eqn.  30). 
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RADIAL  STRESS,  or,  (MPa) 


NUMERICAL  DATA 


This  appendix  contains  a  listing  of  the  experimental  data 
shown  on  the  figures.  All  stresses  (including  the  confining 
pressure)  are  shown  as  negative  when  compressive.  The  data 
sources  are  listed:  references  refer  to  the  reference  list 
of  the  text. 
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WESTERLY  GRANITE 


(a)  Brace  et  aH ,  (1966):  initiation  (dilatation  curves). 

(b)  Brace  et  aj_,  (1966):  failure. 

(c)  Holcomb  and  Costin,  (1986):  initiation  (acoustic  emission). 

(d)  Mogi,  (1966):  failure. 

(e)  Holcomb  and  Costin,  (1986):  intermediate  damage  surface. 


AXIAL  STRESS  (MPa) 

PRESSURE  (MPa) 

REFEREN 

-130.00 

0.00 

-90.00 

0.00 

-120.00 

0.00 

-110.00 

0.00 

-270.00 

-50.00 

-430.00 

-100.00 

(a.) 

-490.00 

-100.00 

-670.00 

-150.00 

-662.00 

-162.00 

-825.00 

-200.00 

-850.00 

0.00 

-218.00 

0.00 

-228.00 

0.00 

-229.00 

0.00 

-680.00 

-50.00 

OO 

-905.00 

-100.00 

-1028.00 

-100.00 

-1210.00 

-150.00 

-1232.00 

-162.00 

-1380.00 

-200.00 

-95.00 

-5.00 

-131.00 

-16.00 

(C) 

-188.00 

-33.00 

-240.00 

0.00 

-390.00 

-10.00 

-490.00 

-17.00 

-590.00 

-27.00 

-780.00 

-50.00 

-920.00 

-77.00 

-1100.00 

-100.00 

-1500.00 

-200.00 

-1900.00 

-310.00 

-212.00 

-5.00 

-268.00 

-15.00 

-343.00 

-33.00 

ce) 

-415.00 

-52.00 

-501.00 

-69.00 
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WESTERLY  GRANITE 


(a)  Brace  et  £]_,  (1966):  failure. 

(b)  Mogi,  (1966):  failure. 

(c)  Schock  and  Heard,  (1974):  failure. 

(d)  Shlmada,  (1981):  Man-nari  granite,  failure. 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


1100.00 

1500.00 

1900.00 

2100.00 


I 


0.00 

-1910.00 

-1590.00 

-6150.00 

-5220.00 

-4390.00 

-3780.00 

-2950.00 


-77. 

-100. 

-200. 

-310. 

-400 


-18 

0. 

00 

166 

227 

9.00 

2.00 

-305. 

-500. 

00 

00 

279 

2.00 

-682. 

00 

•3323.00 

-941. 

00 

•365 

LEE 

-1182. 

00 

3718.00 

-1264. 

00 

3986.00 

-1468. 

00 

4318.00 

-1782. 

00 

4591.00 

-1909. 

00 

4759.00 

-2023. 

00 

0.00 

-500.00 

-250.00 

-3010.00 

-2500.00 

-2010.00 

-1510.00 

-1000.00 


APLITE 


(a)  Brace  etal^,  0966):  Initiation  (from  onset  of  dilatancy). 

( b)  Brace  et  al ,  (1966):  fracture. 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


-260.00 

0.00 

-280.00 

0.00 

-527.00 

-77.00 

-621.00 

-81.00 

-683.00 

-143.00 

-1178.00 

-238.00 

-1283.00 

-283.00 

ZRHSEcftliJfiTiS 

-605.00 

0.00 

-595.00 

0.00 

-797.00 

-77.00 

-831.00 

-81.00 

-1523.00 

-143.00 

-1768.00 

-238.00 

-2263.00 

-283.00 

-2390.00 

-320.00 
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DUNITE 


(a)  Shimada  et  aj_  (1983):  failure  with  much  acoustic  emission. 

(b)  Shimada  e£  ail_  (1983):  failure  with  modest  acoustic  emission. 

(c)  Shimada  et  al[  (1983):  failure  with  little  acoustic  emission. 

(<*•)  Shimada  et  al  (1983):  initiation  (from  onset  of  dilatancy). 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


-190.00 

-724.00 

-1187.00 

-1630.00 

-1970.00 

0.00 

-54.00 

-137.00 

-250.00 

-340.00 

-2375.00 

-435.00 

-2510.00 

(M 

PHMEM.liJliliHH 

PKiiwi 

-4190.00 

-1000.00 

-3900.00 

-1270.00 

-1500.00 

-5460.00 

-1990.00 

Cc) 

-5100.00 

-6250.00 

-2500.00 

-6330.00 

-2750.00 

-40.00 

0.00 

-254.00 

-54.00 

-437.00 

-137.00 

-1135.00 

-435.00 

vP7 

-1970.00 

-770.00 

-2870.00 

-1270.00 
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ECLOGUE 


(a)  Shimada  et  al_  (1983): 

(b)  Shimada  et  al_  (1983): 

(c)  Shimada  e£  al_  (1983): 

(d)  Shimada  et  al  (1983): 


failure  with  much  acoustic  emission, 
failure  with  modest  acoustic  emission, 
failure  with  little  acoustic  emission, 
initiation  (from  onset  of  di latency). 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


-200.00 

-676.00 

-1172.00 

-1600.00 

-1681.00 

-2190.00 

-2373.00 

-2580.00 

0.00 
-56.00 
-172.00 
-240.00 
-251.00 
-350.00 
-453.00 
•  -510.00 

(0) 

-3650.00 

-1020.00 

-4910.00 

-1490.00 

(b) 

-6270.00 

-5830.00 

-2400.00 

-6250.00 

-2490.00 

(0 

-7520.00 

-3010.00 

00 

00 

00 

00 

00 


<f) 


GABBRO 


(a)  Shimada  et  al_  (1983): 

(b)  Shimada  et  al_  (1983): 

(c)  Shimada  et  al  (1983): 

(d)  Shimada  et  al  (1983): 


failure  with  much  acoustic  emission, 
failure  with  modest  acoustic  emission, 
failure  with  little  acoustic  emission, 
initiation  (from  onset  of  dilatancy). 


AXIAL  STRESS  (MPa) 

PRESSURE  (MPa) 

REFERENCE 

-210.00 

0.00 

-383.00 

-53.00 

-769.00 

-149.00 

-240.00 

-997.00 

-247.00 

-1493.00 

-343.00 

-1528.00 

-398.00 

-510.00 

-2210.00 

-760.00 

(b) 

-2840.00 

-990.00 

-3860.00 

-1500.00 

-4210.00 

-1990.00 

(<) 

-5220.00 

■  ll  1  n  l  — 

-303.00 

-53.00 

-697.00 

-247.00 

W 

-943.00 

-343.00 

-1310.00 

-510.00 

-1760.00 

-760.00 

I 

I 
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BUNTE  SANDSTONE 

(a)  Gowd  and  Runnel,  (1980):  failure. 

(b)  Gowd  and  Rummel,  (1980):  initiation  (from  onset  of  di latency). 

AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


SOLENHOFEN  LIMESTONE 


(a)  Heard,  (I960):  failure  (ductile,  axial  stress  dominant), 

(b)  Heard,  (1960):  failure  (transitional,  axial  stress  dominant) 

(c)  Heard,  (1960):  failure  (brittle,  axial  stress  dominant). 

(d)  Heard,  (1960):  failure  (ductile,  radial  stress  dominant). 

(e)  Heard,  (1960):  failure  (brittle,  radial  stress  dominant). 

(f)  Heard,  (1960):  initiation  (from  onset  of  dilatancy). 

AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


-1240.00 

-955.00 

-965.00 

-1010.00 

-667.00 


-500.00 

-300.00 

-300.00 

-300.00 

-150.00 


MARBLE 


(a)  Brace  et  aj_,  (1966):  initiation  (from  onset  of  dilatancy). 

(b)  Brace  et  al ,  (1966):  fracture. 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 

f 

-20.00  0.00 

-115.00  -25.00  (a) 

-179.00 _ -49.00 _ 

-46.00  0.00 

-155.00  -25.00  (b)  j 

-309 . 00  -49.00 


♦ 
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CARRARA  MARBLE 


(a)  Von  Karman,  (1911):  failure  (brittle). 

(b)  Von  Karman,  (1911):  failure  (transitional). 

(c)  Von  Karman,  (1911):  failure  (ductile). 

(d)  Edmond  and  Paterson,  (1972):  failure. 


AXIAL  STRESS  (MPa) 

PRESSURE  (MPa) 

REFERENCE 

-134.00 

0.00 

L«M 

-230.00 

-24.00 

-315.00 

-45.00 

-345.00 

l*) 

-376.00 

-76.00 

-518.00 

-147.00 

-225.00 

tc) 

-735.00 

-295.00 

-356.00 

-100.00 

(*) 

-523.00 

-200.00 

-983.00 

-600.00 

(*) 

-1180.00 

-800.00 

-236.00 

-600.00 

ROCK  SALT 


(a)  Hunsche,  (1984):  failure  (standard  tri axial  test). 

(b)  Hunsche,  (1984):  failure  (cubic  samples,  multiaxial  press) 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


-36.60 

0.00 

-42.00 

-2.50 

-43.00 

-2.90 

-45.90 

-3.00 

-48.70 

-4.00 

-48.50 

-5.20 

-49.50 

-4.80 

—  -50-20 

-4.80 

-52.10 

-3.50 

-54.20 

-5.00 

-58.30 

-5.20 

-58.60 

-6.40 

-64.00 

-8.20 

-68.70 

-10.40 

-71.30 

-10.00 

W 

-69.10 

-10.80 

-76.30 

-15.90 

-20.30 

-81.00 

-19.70 

-81.60 

-19.80 

-82.40 

-19.80 

_ -57.40 

-30.50 

-53.60 

-6.70 

-54.40 

-5.80 

-55.30 

-6.10 

-56. *0 

-6.50 

-60.50 

-9.40 

-64.70 

-10.60 

-65.80 

-11.30 

-66.10 

-10.30 

-68.40 

-10.10 

-64.00 

-8.60 

0>) 

-65.50 

-12.90 

-70.30 

-12.80 

-71.00 

-12.70 

-72.20 

-12.40 

-73.80 

-13.20 

-73.50 

-14.90 

-81.70 

-19.10 

-84.60 

-21 .40 

-86.60 

-20.20 

-88.10 

-18.10 

-96.80 

-25.50 

-98.30 

-29.30 

-99.10 

-37.00 

-108.10 

-41.70 
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ROCK  SALT 

(a)  Handin,  (1953):  failure  (axial  stress  dominant). 

(b)  Handin,  (1953):  failure  (radial  stress  dominant). 


AXIAL  STRESS  (MPa)  PRESSURE  (MPa)  REFERENCE 


-26.30 

-1.43 

-43.00 

-2.23 

-58.04 

-10.72 

-71.27 

-17.41 

-90.93 

-51.16 

-130.77 

-77.78 

-186.72 

-124.86 

-197.70 

-124.79 

-262.29 

-210.40 

-359.91 

-285.23 

-463.24 

-516.67 

-441.91 

-493.32 

-349.63 

-407.64 

-254.03 

-307.98 

0>) 

-240.16 

-249.97 

-187.22 

-276.09 

-147.71 

-202.51 

-69.98 

-99.13 
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